Work, free energy and dissipation in voltage driven single-electron transitions 
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We apply the general procedure presented by Jarzynski in [l[ to the problem of dissipation in a 
volt age- driven single-electron box. We obtain the expression of dissipated work, and find its relation 
to the dissipation Q obtained in 0] . We show that the two quantities are identical in common gate 
protocols where the system makes a transition for sure. 

I. INTRODUCTION 

The non-equilibrium fluctuation relations yield important results for the distribution of dissipation in exper- 
iments where external control parameters are varied non-adiabatically 6-8] . The validity of these relations rests on 
appropriate expressions for work and free-energy. This issue was highlighted in a scholarly manner in [l|. The pre- 
sented general treatment allows one to apply it to a specific system, which in the current case is a single-electron box 
(SEB) [3, [13], see Fig. QJa). Single-electron transitions [ll[ yield a suitable test-bench of fluctuation relations, where 
reliable statistics can be collected under stable experimental conditions [l2j • In a SEB, a tunnel contact, schematically 
a split box on the left in Fig. QJa) admits electrons to enter or leave the island in the middle at rates determined by 
standard tunneling expressions. The characteristic electrostatic energy of the box can be made large (small size, low 
temperature) such that the number of electrons on it is limited to only two neighbouring values in a given control 
parameter (gate) interval. Our analysis is, however, not restricted to such two-state dynamics. 

The quantity Q introduced in Ref. Q is the dissipated energy in the tunneling transitions only, and thus it represents 
the heat generated. In the present paper we demonstrate mathematically that the dissipated work, W — AF, relevant 
in fluctuation relations [3, 0, [j| , consists of the sum of Q and the electrostatic energy stored in the circuit during the 
driven evolution. 
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II. WORK AND FREE-ENERGY IN A SINGLE-ELECTRON BOX 

We discuss the hamiltonian H of the SEB, yielding the work W in a given protocol of the external drive voltage, 
and the equilibrium free-energy F at a given voltage bias, along the general prescription given in [l|. 

The "bare" hamiltonian Hq of the SEB is given by the bare electrostatic energy of the capacitors in the circuit, 
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Ql + Ql 



a) 



where Cj is the capacitance of the tunnel junction, and C g is that of the gate, see Fig. [DJa); Qj is the charge at the 
junction and Q g is that of the gate capacitor. 

The generalized force of the SEB is the gate voltage V g , and the generalized coordinate is its charge Q g . With the 
help of these, we can express the transformed hamiltonian H as 
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FIG. 1: (a) The single-electron box with a voltage source for controlling the gate charge, (b) The voltage source represented 
by a large capacitor. 
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We give a proper justification of this basic transformation in Section IIIIl We next write down H as a function of the 
integer number n of electrons on the island and the normalized gate voltage n g = —C g V g /e. The gate charge is given 
by Q g = C g (V g — (p a )- Here <p a is the island potential 

4>a = -jj-(n-n g ), (3) 
and Cs = Cj + C g is the total capacitance of the SEB. By elementary analysis we then have 

H(n,n g ) = E c n 2 - 2E c nn g - ^-(Cg 1 - C^)n% (4) 



where Ec — e 2 /(2Cs) is the charging energy of the SEB. 

The work performed in a driven process can be written in general as W = J Q^-\dt, where A is the control parameter 
and the dot refers to the derivative with respect to time t Thus the work for our system reads 



/ ^ g dV g , (5) 



where the integral is extended over the gate excursion between the end points. Direct differentiation of Eq. (j4]) yields 

ML = —O 

dv„ ^9 



— ~Q g and the proper "thermodynamic" work is 



W = -jQ g dV g (6) 

in full analogy to the result in The "standard" expression of work is, on the other hand, the integral of the force 
along the position, i.e., 

W = jVgdQg. (7) 

As discussed in Ref. [lj, the two expressions ([6]) and (JT)) are identical for a cyclic trajectory where the perturbation 
V g is turned on and then off, as can be verified by elementary partial integration. 
With the notations we have used above, we find 

W - -2E C / ndn g - -(CZ 1 - C^){nl B - nl A ), (8) 



where n g J \ and n g B are the gate positions in the beginning and at the end, respectively. 

Next we find the change in free energy, again according to the standard procedure [lj. For completeness, we allow 
the system to occupy any integer valued charge state n although it is exponentially unlikely for the system to visit a 
state for which \n — n g \ 3> 1 / V p Ec , where /3 is the inverse temperature of the bath. For writing down the partition 
function Z(n g ), we first regroup the Hamiltonian as H — Ec{n~n g ) 2 — e 2 / (2C g )n 2 , and then perform the summation 
over possible n states as 

Z(n g ) = e' fjH ^ = exp (^-) £ e ^ c (n-n a )^ (9) 

n y n 

The remaining infinite summation defines a special function which we will refer to as RpE c { n g)- The function is 
periodic in n g with a period of 1, which is a manifestation of the fact that there is no absolute reference point for the 
charge number n. 

The change of free energy for the given gate sweep n g ,A ~> n g,B is then 

AF = F (n gtB )-F(n g , A ) = -r 1 hi^f4- 

= -^KB-<,)-r 1 lnf^l. (10 ) 

The dissipated work, W — A_F, the central quantity in non-equilibrium fluctuation relations, can then be written 
combining Eqs. © and (fTU)) as 

W~AF = 2E C (n g - n) dn g + /T 1 h ^^j . (11) 

Jn aA ttj3EA n g,A) 
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In the basic gate protocol with n g .A = 0, n g .B = 1, used in the recent experiment by Saira et al. [12j . the expression 
for dissipated work assumes the form 

W- AF = E c (l - 2 / ndn g ). (12) 
Jo 

The same expression is obtained also if only states n — and n — 1 are allowed in the calculation of the partition 
function. This is the appropriate choice for, e.g., experiments in the high charging energy regime fiEc 3> 1. 



III. FIRST-PRINCIPLES DERIVATION OF THE CLASSICAL HAMILTONIAN 

We give an elementary derivation of Eq. (j4|) based on Ref. [T3|. Following the general procedure, the gate voltage 
source in Fig. QJa) is represented by a large capacitance Cs in which is stored intially a large charge Q s such that 
Q%/Cg — V g . First, one writes the Lagrangian in terms of the node fluxes (t> a .b '■— f_ 4>a,b{t')dt' , where <p a ^ refers 
to the potential of nodes a, b as defined in Fig. [IJb). The result is 

£fuu - + \c g (J> a - 4> b ) 2 + has fa. (13) 

It is important to distinguish between <pb, referring to the variable voltage of node b, and V g , the initial voltage of the 
source capacitor Cs- The generalized momenta are obtained with the usual relation q a ,b — g^ fu " , and are equal to 
the sum of charges at the nodes a, b, 



Qa = Cj<j> a + C g (<j> a - fo) (14) 
q b = C g {4> b - fa) + C s j>b- (15) 



The Hamiltonian Hf u \\ is obtained as 



H{ u n = (j) a q a + 4>bQb - £full- (16) 

To obtain a reduced Hamiltonian describing the dynamics of the box charge exactly, one needs to subtract the energy 
of the capacitor Cs from the full Hamiltonian. Subtracting the instantaneous energy ^Cs^, one obtains 

fffuii - hosft = hoAl + \c g {k - « 2 , (17) 

which equals Hq of Eq. ([1}, as it should. On the other hand, subtracting the energy hCsV g , which was stored initially 
in the capacitor should yield the Hamiltonian H of Eq. @ in the limit Cs — > oo [13| . We obtain 
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Hfiui - i;C s V g 2 = H + -C s (<ft - V?). (18) 



We next analyze the difference ^Cs{<Pl — Vg). The charge is conserved on the node b. Thus the gate charge at an 

arbitrary time (assuming it is zero initially) reads Q g = Cs(V g — 4>b). Inserting this charge conservation into the said 
difference yields 

±C s {j>l-V g 2 ) = ~Q g {j> b + V g ). (19) 
In the limit Cs — > oo, we have 4>b — > V g , yielding 

H - H = ho s {fa - V g 2 ) = -Q g V g (20) 



in accordance with Eq. (|2]| . 
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FIG. 2: Graphical interpretation of Eqs. (|23[) - 126} to compare W — AF and Q along gate trajectories n g : — > 1 with different 
outcomes. Outcome A, n : — > 1, is the expected (and in the recent experiments [lj] the only) choice, where the system makes 
the expected transition. There W — AF = Q. In the unlikely outcomes of B, C, and D, the electron numbers in the beginning 
and at the end are (0, 0), (1, 1), and (1, 0), respectively. The dashed lines in the middle panel give full simulated Q distributions 
for different types of outcomes, where the capital letter refers to the four outcomes A - D, weighted with their probability. The 
simulations were made for a / = 5 Hz sinusoidal drive using the parameters of Ref. Ec/ks = 1.94 K, temperature of the 
experiment T = (fcs/?)" 1 = 214 mK and the junction between superconducting (aluminium, A = 218 ^teV) and normal metal 
leads had a resistance Rt = 100 Mfi. The full distribution of W — AF, the solid line, is indistinguishable from that of A of 
the dominating outcome. For the full distribution the horizontal axis is naturally replaced by (W — AF)/Ec- 



IV. COMPARISON OF THE DISSIPATED WORK W — AF AND THE DISSIPATION Q IN TUNNELING 

The dissipated energy in tunneling transitions along the gate drive was obtained in Ref. [2| as 

Q = E c ^2±(2n g ,i-l), (21) 

i 

where the sum is over all back and forth tunneling events i along the ramp, and ± refers to the direction of the jump: 
+ into the box, and — out of the box. 

We now compare the results (Tj"2"|) and (j2"Tj) for different outcomes of the experiment. For convenience, we integrate 
Eq. (fl~2"|) by parts and obtain 



W-AF = E C (l-2n f + 2j2±n g ^ = E c (1 - m - »/) + Q. 



(22) 



where is the charge state at the beginning (end) of the protocol and the sign in the summation is chosen as in 
Eq. (fT2")l . For the most obvious case, i.e., the "successful" trajectories n : — > 1 (meaning that the charge tunnels in 
the ramp as expected, perhaps with several intermediate back and forth transitions), we then have 



(W - AF) ^ =E C J2 ±( 2n 9.« - 1) = 0- 

i 

For trajectories n : — > (even number of jumps) 

(W - AF) ^o = E C [Y^ ±(2n 9il -1) + 1] = Q + E C 

i 

For the other types of trajectories we have 

(W - AF)^ = E C ^2 ±(2%^ - !) = Q 



and 



(W - AF)i_n = E C [J2 ±( 2 %,i -1)-1] = Q-E C 



(23) 



(24) 



(25) 



(26) 
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V. INTERPRETATION 

Figure [2] gives a graphical interpretation of the comparison in Eqs. (f23]) - (|26p . The two quantities, the dissipated 
work, W—AF, and dissipation Q are directly related. The difference between them is that Q gives the pure dissipation 
in the tunneling events, which is then released as heat typically to the phonon system, whereas W — AF is the sum 
of Q and the energy stored in the electronic system in the gate ramp. In the properly designed experiment where the 
transition takes place in all realizations (fiEc 3> 1) [l^, the two quantities are equal. 
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